Density waves and jamming transition in cellular automaton models for traffic flow 
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In this paper computer simulation results of higher order density correlation for cellular au- 
tomaton models of traffic flow are presented. The examinations show the jamming transition as 
a function of both the density and the magnitude of noise and allow to calculate the velocity of 
upstream moving jams. This velocity is independent of the density and decreases with growing 
noise. The point of maximum flow in the fundamental diagram determines its value. For that it is 
not necessary to define explicitly jams in the language of the selected model, but only based upon 
the well defined characteristic density profiles along the line. 



I. INTRODUCTION 

Recently, the examination and modeling of vehicular traffic has become an important subject of research - see Q- 
H and references therein for a brief review. In the microscopic approach to the traffic flow problem, the cellular 
automaton introduced in |6j reproduces important entities of real traffic, like the flow-density relation or stop-and-go 
waves. Beside the realization of some basic requirements to such a model it can be efficiently used in computational 



investigations and applications j?]]- |10|. Fundamental analytical and numerical examinations enclose exact solutions 
for certain limits and mean- field approximations fllTll , the jamming transition fll2||- JT7[ or the effects of perturbations 
and the occurrence of meta-stable states fisf| - [pOfl , e.g. We investigate the density waves and the separation in 
free-flow and dense regions by means of the density-autocorrelation function. It enables us to trace back the spatio- 
temporal evolution of jams which are stable during the measurement time, on condition that jams emerge. It should 
be noted that the probability for a jam to survive decreases with the simulation time pl[ | in a system without a 
clear phase separation between congestion and free-flow. But the duration of a simulation is sufficiently shorter than 
these time periods. By this method it is superfluous to give an explicit definition of what a jam is and which cars 
are belonging to the jam. Therefore this method can be theoretically used for every traffic flow model where density 
profiles are available. As an example we apply this method to cellular automaton models. In this context we will 
report and discuss several aspects of the underlying model and their slow-to-start modifications (Section ||): The 
jamming transition shows up by varying both the global density p and the global noise p; the jam velocity can be 
derived dire ctly from the density-autocorrelation function and is closely related to the global flow-density relation 



(Section III). It is not the goal of the paper to discuss the jamming transition with regard to the criticality or the 



sharpness of this issue, but we apply the method of correlation function from a more practical point of view. 



II. THE MODEL 



Within the framework of this paper we only consider a one-dimensional ring of cells. The cells are either vacant or 
occupied by a vehicle labeled i. Its position is Xi and its discrete velocity is Vi G [0,u ma3; ]. The gap denotes the 
number of empty sites to its leading vehicle. The rules for a parallel update are 

• Acceleration with regard to the vehicle ahead: v[ <— min(«i + 1, gi,v max ), 

• Noise: with a probability p do v'( <— max(uj — 1,0), 

• Movement: Xi *— xi + v'/ . 

The investigated systems consist of L cells and N vehicles, the global density is p = N/L. The flow is defined as 
J = (v)p with the mean velocity (v) = ^Vi/N . In the following the Nagel-Schreckenberg cellular automaton model 
H defined through the above set of rules is denoted by SCA. 

We extend our studies on further modifications of the SCA, namely on models with slow-to-start rules. For the 
model with velocity-depending randomization (VDR) [jl9],^0) we set p(vi — 0) = Mm(p + pvdr, 1)- That leads to a 
reduced outflow from a jam. Note that Vi is the velocity before the first update step is performed. The other modified 



model under consideration is the T 2 -model introduced by Takayasu and Takayasu 22 1. Here the headway gi of a 



1 



vehicle i controls the acceleration: standing vehicles with a headway gi = 1 only speed up with a probability 1 — p 
with p = Min(p+pT2, 1), whereas for all other ones the rules are unchanged. Unlike the SCA with similar parameters, 
both models exhibit a different behavior in the vicinity of the point of maximum flow (p m g X = p(Jmax), J max)- They 
are capable to generate metastable states in the adiabatic approach (for details see (£0|), i.e. one finds two branches 
of J(p) in a small density interval. Additionally, for sufficient small p one can found a clear separation of of dense and 
free-flow regions in a space-time plot. We summarize the effective deceleration probabilities of the applied models: 

const. 

( Mm(p + pvdr, 1) Vi = 

1 p otherwise (1) 

| Min(p + p T 2,l) Vl = 0A gi = l 
\ p otherwise 

Actually, other definitions of p are conceivable, but, for our purpose, we decided to use only the above notations. 
Primarily, it was done for modeling moving vehicles with similar properties and to scan the parameter space by 
varying only p. Both pvdr and px? are of any value but fix. 



SCA: p = 
VDR: p = 

T 2 : p = 



III. SIMULATION RESULTS AND THEIR DISCUSSION 

The density waves are moving upstream and can be easily observed in a space-time plot PJ23]| , one finds separation 
of dense and free-flow regions. For the measurements it is necessary to introduce the mean local density pi(k,t) of 
the cell k at time t: 



ft(M) = rX>*+i(f) (2) 



i=0 



with 

Vk+i{t) = 



1 if site k + i is occupied at time t 
otherwise 



The parameter A denotes the length of the interval on which the local density has to be computed, ft should satisfy 
the condition Ao C A < L with a characteristic length scale Ao- For the determination of the jam velocity Vj we 
use the generalized T-point-autocorrelation function of the density 



T-l 



C v; (r = VjVAT, t) = (JJ Pi{x + VjrAT, t + rAT)) L (3) 

T = 

with the supposed jam velocity Vj G [—1,0]. By varying Vj one finds a largest Cy*(r,r) (Fig. [I]). AT is the time 
interval between two single measurements which contributes to (|3|). Sufficiently large values of AT are necessary to 
observe a macroscopic motion and to determine Vj with an adequate accuracy. Unless otherwise mentioned, we set 
L = 10 4 and AT = 10 2 in order to exclude any finite size effects. Usually, we are averaging over 20 simulation runs 
with v m ax = 5 and A = 30. 

As pointed out in Fig. [I] one has to adjust thoroughly the parameters T and AT. If T is of the order of magnitude 
of AT then the uncertainty of the measurement covers the interesting signal and Cvj vanishes. This problem becomes 
more serious while approaching p* . In this region the calculations are complicated additionally due to large fluctuations 
of Cv.j itself. 

The jam velocity Vj depends on the deceleration probability p (Fig. ||) and is related to global quantities as it will 
be shown later. We checked this for a variety of parameters, but could not notice any remarkable deviations among 
the diverse sets of data. The modifications of the SCA, VDR and T 2 , yield a different behavior. The absolute values 
of the jam velocity in the models VDR and T 2 are smaller than that of the SCA implying that both modifications 
are characterized by a lowered outflow from a jam that, in turn, reduces the jam velocity. 

The results of simulations using the VDR are basically shifted towards smaller |Vj|'s. pvdr can be recognized on 
the ordinate at p = and on the abscissa at Vj = 0, for large p's a total deadlock occurs, i.e. it is highly unlikely or 
even impossible that a stopped car speeds up again. In order to find out how Vj is related to p we investigate the 
mean waiting time t w for p < 1 expressed through an infinite series: 



2 




FIG. 1. The peak of the density-autocorrelation function enables to estimate the jam velocity Vj. The standard deviation 
of the symmetrically assumed CV ; is depicted in the inset. For large ratios T/AT the variance shrinks a lot, therefore it might 
happen that the autocorrelation function vanishes and inhibits the estimation of depending quantities. In the vicinity of p* 
this sensitivity is more pronounced (p = 0.4). 
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FIG. 2. The jam velocity as a function of p, the error bars are within the symbol size (p = 0.4). 
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FIG. 3. Vj vs. p for p — 0.5. Beyond p* , especially for 0.2 < p < 0.4, Vj(p) can be assumed to be constant. The inset reveals 
the large fluctuations of Vj nearby p* which are up to the order of magnitude of Vj itself. 

oo 

t w = l(l-p) + 2(l-p)p + 3(l-p)f ■■■ = (!- p) Vnf-^- = 

^ 1 ~P 

=>\Vj\ = ^- = l-p G [0,1 -Pvdr] (4) 

and is exact for p = 0. For large values of pvdr one obtains a good agreement, whereas for small pvdr the jam 
velocity is overestimated. This is due to the so-called sub-jams which emerge downstream a wider jams and cause a 
reduction of \Vj\. 

The results drawn from simulations using the T 2 model show no deadlock situation for any p < 1. Starting with 
p = 1, one can hardly distinguish between the simulation results of SCA and T 2 . Especially, this is valid as long as 
p T -i > 1 — p. Similar to (^) one can estimate 

t w = 1(1 -p) + 2(1 - p)p + 3(1 - p)pp + 4(1 - p)pp 2 ■ ■ ■ 

n=0 ^ 

Note that for small p T i -values Vj is overestimated for all values of p, which, in turn, can be traced back to the 
occurrence of sub-jams. With increasing p T i even for small p it is required to set p = 1. Again, it is \Vj\ — and 
two special cases can be described by 

\Vj\(p = 0) = —±— and \Vj\( PT 2 - 1) = . (6) 

Obviously, the measurements reveals several density regimes. Below p* the vehicles move independently, i.e. there 
are no correlations between them. For p > p* upstream moving density waves can be detected by means of In 
the vicinity of p* the jam velocity reveals large fluctuations (Fig. |3|), which are due to the recurrent emergence and 
dissipation of jams. But beyond p* Vj is nearly constant. Within the interval where density waves are to be expected 
it is obvious that Vj is independent of p, since the outflow from a jam is independent from the global density. 

So far, we applied the autocorrelation function (^) to determine the jam velocity. But this quantity itself indicates 
the two different phases separated by noise p* or density p* (Fig. [I| and ||). Varying p leads to a transition while 
crossing p* . Its clarity strongly depends on T/AT, for insufficient ratios a plateau at Cvj(p) occurs. To elucidate it 
we used a modified autocorrelation 

T-l 

Cvj(r,r) = ((H Pl (x + Vj*rAT,t + tAT))V t ) l . (7) 

T = 

The other transition takes place while crossing the density p* (Fig. ^1). For p < p* one finds empty regions on the 
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FIG. 4. The transition from free flow to congested flow can also be obtained in the behavior of the modified autocorrelation 
function Cy } (fjj) by varying p (p = 0.073). The transition is smeared out due to finite size effects and systematic errors in the 
determination of CVj- 




0.0 0.2 0.4 0.6 0.8 1.0 

Global density p [veh/site] 



FIG. 5. Plot of the the autocorrelation function C\rj {p) \v m ax ■ Below a signified density the autocorrelation function vanishes 
due to the absence of jams. The inset zooms into the region p ~ p* for the modified autocorrelation Cvj (R) (p = 0.5). 
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FIG. 6. Plot of the the autocorrelation function CV ; (v ma x)\ P (p = 0.5). 

road of the order of magnitude of A, and therefore Cvj completely vanishes. On the other hand, for p > p* stable 
congestion emerge. The same jam can be detected at U as well as at tf = U + tAT located at x{ti) — \ Vj\tf. In this 
context, p* can be denoted as the density, at which stable jams emerge and separates the density regime as it is to 
be seen in the inset of Fig. [|. For a fixed density and a varying v max (Fig. [jj) the relationship can be estimated as 
Cv.j {vmax) °c p- Nevertheless, the quality of these data does not allow a correct classification of the transition between 
the free-flow and the dense region. Above all, the sensitive dependencies on the many adjustable parameters of this 
method seem to prevent a more accurate consideration of the interesting interval of density. 

How is the jam velocity related to other macroscopic quantities? In the steady state the dynamics are characterized 
by an equilibrium of out-flowing vehicles and vehicles attaching the jam from behind. The more frequent vehicles 
join the jam the faster the jam moves upstream. If we neglect any effects due to meta-stability then the free-flow 
region can be assumed to be located in the vicinity of the point of maximum flow (p ma x, Jmax)' Hence, the velocity 
of attaching vehicles is (v a tt) — Jmax/ Pmax- The mean distance between the upstream tail of the jam and the next 
vehicle is (g) = p m \ x — 1 and the temporal distance therefore reads 

M att = M- «. Vj ^-^-<0. (8) 

\Vatt) Pmax — 1 

This is confirmed by the simulation results depicted in Fig. It means that Vj is determined by the slope of the 
congested branch (p > p m ax) in the fundamental diagram. This can also be verified for the modifications VDR and T 2 
(Fig. ^J). The small deviations from the data rest upon a difference between the outflow of the jam and the maximum 
global flow in the considered systems, but also in the above made assumption of the equilibrium. Actually, the lowered 
outflow from jams observed in the models VDR and T 2 in comparison to the SCA is also reflected by (Js|) . 

Concluding, this knowledge enables a calibration of the SCA. Besides the approach of the fundamental diagram 
derived from empirical data a further point of interest is the velocity of upstream moving jams (« —15 km/h on 
German highways p4[) - but according to (|^) all information is accumulated in the fundamental diagram, namely in 
the second characteristic slope of J(p). For the SCA one can set v max = 5 and p — 0.2 • • ■ 0.3 to adapt the simulation 
to this empirical jam velocity. 

IV. SUMMARY 

We investigated the cellular automaton model for vehicular traffic in order to get information about the density 
waves and their velocity. Beside the standard Nagel-Schreckenberg cellular automata we also included two slow-to- 
start modifications (VDR and T 2 ). Both resemble the SCA except the rules for standing vehicles. Loosely spoken, they 
result in a lower flow downstream a jam and a clear phase separation for certain density regimes. For the determination 
of the jam velocity we used the density-autocorrelation function Cvj(t,t). Despite the high computational efforts 
(0(L 2 )) this method was suitable to be applied for our calculations. Moreover, a definition of jams is not necessary 
and therefore the method can be applied to every model that provides density profiles along the road. 
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FIG. 7. The jam velocity can be explained by (H): Vj = J ma x / (Pmax — 1) 
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FIG. 8. Vj (filled symbols) as well as Jmax / [pmax — 1) (opaque symbols) are depicted for all used models. The small 
deviations are related to the discrepancies between the outflow of a jam and the global maximum flow. 
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The quantity C'vj (r, t) reflects the two different phases and depends on the global density p. The density regime 
is separated by p* . For p < p* no jams can be detected by the applied method, whereas for larger p the system 
is dominated by sequences of dense and free- flow regions, where Cvj remains finite and permits to estimate p* . At 
this point a transition to the congested region takes place. Both the local length A and the number of calculations 
T have large influence on Cvj- Further statements regarding the transition cannot be given due to the numerical 
insufficiencies and accuracy. 

The jam velocity can be derived directly from CV,(r, r). For sufficiently large p the absolute value of Vj is a 
continuous and descending function of p, but depends on the considered model. The differences between the models, 
especially for p — > and p — ► 1, could be explained by waiting time arguments. The jam velocity is essentially 
expressed through p max and J max , irrespective of the model considered here. 
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